Dependence of nuclear spin singlet lifetimes on RF spin-locking power 
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Abstract 



We measure the lifetime of long-lived nuclear spin singlet states as a function of the strength of the RF spin-locking 
field and present a simple theoretical model that agrees well with our measurements, including the low-RF-power 
regime. We also measure the lifetime of a long-lived coherence between singlet and triplet states that does not require 
a spin-locking field for preservation. Our results indicate that for many molecules, singlet states can be created using 
weak RF spin-locking fields: more than two orders of magnitude lower RF power than in previous studies. Our 
findings suggest that in many biomolecules, singlets and related states with enhanced lifetimes might be achievable 
in vivo with safe levels of RF power. 



1. Introduction 



The spin-lattice relaxation time T\ is a limiting factor for a broad class of NMR experiments in which nuclear spin 
polarization or order needs to be preserved or transported . However, long-lived nuclear spin singlet states with 
lifetimes up to 37 T\ have recently been measured in thermally-polarized samples [5H10j . These long-lived states 
have been used to study slow processes such as diffusion, chemical exchange, and conformational dynamics in vitro 

\niuM- 

Nuclear spin singlet states are typically created from pairs of coupled spin-1/2 nuclei with equal or near-equal 
resonance frequencies. The most well-known example is the H2 molecule, which exists in two forms: para-H2, with 
the singlet eigenstatc \S ) = (| f-i) - I 4-t))/V^i and ortho-H 2 , with triplet eigenstates |T_) = | tt)> \ T o) = (| tl 
) + I ~i-t))/v2 ; and \T + ) = \ H) 14J. In our notation f represents a spin aligned with the applied magnetic field, B , 
while J. represents a spin that is anti-aligned. In the singlet state, the net nuclear spin is zero, and there is no net 
magnetic dipole moment. Hence interactions with the environment are weak and the rate of interconversion between 
singlet and triplet states is very slow, often much slower than the spin-lattice relaxation rate 1/Ti. On the other 
hand, the triplet states have non-zero magnetic moments and couple strongly with the environment. Relaxation 
among the triplet states occurs on the timescale T±. 

Recently, Levitt and colleagues demonstrated a general technique for the creation of singlet states for pairs of 
magnetically inequivalent nuclei [7J [T3] . Such inequivalent nuclear spins cannot form ideal, long-lived singlet states 
as outlined above: the different local environments of the nuclear spins leads to rapid conversion to the triplet 
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state and thus coupling to the environment. Nonetheless, as Levitt et al. showed, a properly designed RF pulse 
sequence can prepare a singlet state, which is then preserved from triplet interconversion by the application of a 
continuous resonant RF field. This "spin-locking" field forces the average Hamiltonian of the two nuclear spins to 
be effectively equivalent. While this technique is applicable for a large variety of molecules, the large continuous 
RF power employed in spin-locking experiments to date implies an RF specific absorption rate (SAR) that is likely 
prohibitive for animal and human studies |16[ I17j . 

In this paper, we report measurements of singlet state lifetime for a variety of organic molecules and as a function 
of RF spin-locking field strength. We find that the measured RF power required to preserve a singlet state agrees 
well with the predictions of a simple theoretical model with inputs from the molecule's NMR spectrum. We also 
present measurements of a singlet-triplet coherence with an extended lifetime that does not require the use of RF spin 
locking for preservation. Moreover, our findings demonstrate that for many molecules of interest, singlet lifetimes 
many times longer than T\ can be achieved with much weaker RF spin-locking fields than have been used to date - 
more than an order of magnitude smaller than in previous studies - leading to both an RF power and an SAR more 
than 100 times lower. This result suggests that in vivo application of long-lived singlet NMR might be possible in 
biomolecules with the appropriate properties, despite limitations imposed by RF SAR. 



2. Theory 

2.1. Singlet Relaxation Mechanisms 

Many relaxation mechanisms are forbidden by symmetry from converting the singlet state to triplet states. For 
example, magnetic dipole-dipole interactions between the singlet's two spins cannot couple the antisymmetric singlet 
state to the symmetric triplet states. Since this intra-pair interaction is often the dominant driver of relaxation, 
the typical result is a singlet population with a lifetime T$ many times longer than the spin-lattice relaxation time 
T\. Dipole-dipole interactions between the singlet spin pair and more distant spins can also lead to relaxation, but 
the singlet is protected from dipolar fluctuations common to both singlet spins: i.e., in the far-field the net dipole 
moment of the singlet is zero. Thus singlet-state relaxation must instead occur through differential interactions on 
each spin of the singlet; these include chemical shift anisotropy (CSA), spin rotation (SR) due to collisions between 
molecules, and magnetic quadrupole interactions with a third spin |18H21j . 

Since these three singlet relaxation mechanisms respond differently to the applied static magnetic field and temper- 
ature, the dominant effect can be determined experimentally. CSA has a strong magnetic-field dependence with a 
lifetime scaling as Tqsa oc 1/Sq 18J. Spin rotation collisions result in a lifetime scaling non-linearly with temperature 
a s Tsr oc exp(E/kBT) oc -j^p [18]. Magnetic quadrupole interactions, on the other hand, produce a singlet lifetime 
scaling linearly with temperature in the extreme-narrowing regime (when molecular rotation rates are much greater 
than the Larmor frequency). As demonstrated below, for the molecules used in the present study, the quadrupolar 
mechanism dominates singlet-state relaxation, i.e., T$ ~ Tq. 

Magnetic quadrupole relaxation results from the two spins of the singlet interacting differently with a third spin. 
This relaxation mechanism was modeled at high magnetic field by Tayler et al. [21 j . who derived an expression for 
the enhancement of the singlet lifetime^ 

Ts = 3fe?2 

Ti 2 J2(bi 3 + b% - hj b 2j (3 cos 2 0„ a - 1)) ' 

J">2 



*Ref. 21 contains a typesetting error in which the summation has been taken over the whole expression rather than only the 
denominator. 



2 



a Two weakly-coupled spins b Strong scalar coupling c Strong resonant 




Figure 1: A system containing two spins with resonant transition frequencies ui and V2 can be represented by four spin-pair eigenstatcs. 
(A) For weak spin coupling, relaxation occurs via magnetic dipole-dipole interactions due to zero-, single-, and double-quantum transitions 
(with rates Wo, Wi, and W2). (B) If the two spins are coupled strongly by scalar coupling (J >> Av, the chemical shift splitting), the 
spin-pair eigenstates are mixed into singlet and triplet states. The triplet states still interact via dipole-dipole interactions (solid lines), 
but singlet-triplet transitions are forbidden (dashed lines). (C) Mixing of the spin-pair eigenstates can also be induced by a strong RF 
spin-locking field (v n > 5Af), which produces a singlet state and three mixtures of triplet states. Singlet-triplet transitions are again 
forbidden. The new triplet states exhibit new transition rates W[ and Wj. 

Here spins 1 and 2 compose the singlet while j represents another nearby spin; bjk = 7 2 /r| fc is a measure of the 
dipolar coupling strength between spins; and 4>ij% is the angle between the vectors connecting 1 with j and 2 with 
j. In principle, there is no limit to the singlet lifetime enhancement given the proper molecular geometry. However, 
in practice other relaxation mechanisms gain importance if magnetic quadrupole relaxation is highly suppressed. In 
a previous study, Equation 1 was found to agree well with measurements of singlet-state lifetimes using high RF 
spin- locking power [21] . 



2.2. RF-Power Dependence 

A detailed theoretical analysis of the singlet's lifetime during RF irradiation has been given by Pileio and Levitt, 
who performed exact numerical calculations for the relationship between singlet lifetime and RF power [18] . Here, 
we develop an approximate model that leads to a simple calculation of the measured singlet lifetime at a given 
spin-locking frequency. Our model can easily be fit to measurements of singlet lifetime at a number of RF field 
strengths so that the maximum singlet lifetime can be extracted. 

A pair of spin-1/2 nuclei creates a system with four spin-pair eigenstates as shown in Fig. 1A: | "ft), | f4-}> I 4-t); ar± d 
| H). In general, there is mixing between the two eigenstates with no net z-component of spin (| fj,) and | 4-t)) if the 
spin states are degenerate, or if scalar coupling between the spins is strong compared with any difference between 
the individual spin transition frequencies {y\ and v<x). As a result of such mixing, the spin-pair system is described 
by one sing let state, |5 ) = (| U) ~ I it))/V2, and three triplet states |T_) = | ft), \T ) = (| t4-> + I lt))/V2, and 
\T + ) = | 44-); with the energy levels shown in Fig. IB. In this case, spin polarization cannot be transfered to the 
singlet state from the triplet states via an RF pulse sequence because Hrf\Sq) = 0. 

For many molecules of interest, however, chemical shifts induce a difference between spin transition frequencies 
(A^ = \i>i — 1/2D that is much larger than the scalar coupling, and hence there is little mixing of the bare spin- 
pair eigenstates. In this case unitary transformations, via the RF pulse sequence shown in Fig. 2A, can transfer 
initial thermal spin polarization to the singlet state with at most 50% efficiency by creating the singlet-enhanced 
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Figure 2: (A) The singlet-enhanced superposition state pst is initialized via a 3-pulse preparation sequence. (B) A long-lived coherence 
(LLC) between the singlet and triplet states Pllc> 3 s weu 33 mixtures of the singlet and LLC states, are initialized via a 2-pulsc 
preparation sequence. Both preparation sequences are followed by application of the RF spin-locking field and a signal acquisition pulse. 



superposition state [15] 



Pst =\T )(T \-\So)(So\ 

=| tout I + UtXtt 



(2) 
(3) 



Similarly, the RF pulse sequence shown in Fig. 2B can transfer initial thermal spin polarization into a long-lived 
coherence between the singlet and triplet states with density matrix 



Pllc =\S Q )(T \ + \T Q )(S \ 

=mxm-ut)ut 



(4) 
(5) 



A similar long-lived coherence has been previously studied as a way to extend T 2 



For both pst and pllCi population in the singlet state will be rapidly interconverted with the triplet states on a 
timescale ~ 1/Av, providing strong coupling to the environment and thus rapid relaxation to the thermal state. 
However, a strong on-resonance RF field has been shown to be effective for maintaining the singlet-enhanced super- 
position state pst [3 US] . When the spin-locking RF field is set to the average resonant transition frequency of the 
two spins, the Hamiltonian in the bare spin-pair basis becomes 



n = h 



2 



2 



2 


2 







V1+V2 

2 



(6) 



Here v n is the effective spin nutation frequency due to the RF field, which drives single-quantum spin transitions. 
Also, we assume scalar coupling is weak and neglect its contributions to the Hamiltonian. 
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Diagonalizing this Hamiltonian yields four spin-locked eigenstates given by 



= * sin0(| tt> + I U)) + I cos 0(1 U) - | It)) 

+~(it;> + ut» 
\<t>o)=^(\n)-\u)) 



V2 



cos0(|tt) + |W» 



+_sinff(|t;>-Ut)) 



Isine(|tt> + |U»-^cose(|n)-Ut» 



1 



+2(in> + ut». 

The mixing angle is controlled by the ratio of the spin nutation frequency to the chemical shift splitting: 

2v n 



9 = arctan 



Av 



(7) 

(8) 

(9) 
(10) 
(11) 



At very large nutation rates (y n » Az/), i.e., high RF spin- locking power, the spin- locked eigenstates simplify to 



i0+> =^(in> + ut) + itt> + iu» 

=-L|T ) + i(|T_) + |T + )) 
\<h)=j=Qtf)-\U)) = j=(\T-)-\T+)) 

\M=^=(\n)-\m = \s ) 
\<i>-)=U\n) + \it)-\n)-\u)) 



~V2 



|T )--(|T_) + |T + )). 



(12) 
(13) 
(14) 

(15) 



Note that the spin-locked singlet state \<f)g) corresponds to \Sq) in this limit of large spin nutation (i.e., large RF 
spin-locking field), whereas the three spin-locked triplet states are each mixtures of eigenstates \Tq), \T + ) , and |T_). 
In this case, the initial state psr is described well by Equation 2. The singlet, |£o), is well-protected by the RF 
spin- locking field, and after a short initial period during which the triplet states equilibrate, the remaining \Sq) 
component relaxes exponentially with the characteristic time predicted by Equation 1. 

In the high-RF-power regime, the long-lived coherence pllc is a sum of coherences containing \<p+), \4>-), and \<f>s), 
which experience decoherence due to both dipole-dipole interactions and inhomogeneities in the RF spin-locking 
field. 

If instead very small RF spin-locking power is applied (y n « Al>), the singlet component of psr rapidly interconverts 
with the central triplet state, \Tq). When no RF power is applied, psr is a zero-quantum coherence that precesses in 
the xy-plane, with a lifetime up to 3.25 T\ if inter-pair dipole-dipole interactions are the sole relaxation mechanism 
[2~3] . The addition of a small amount of RF power quickly decreases the lifetime of the psr coherence because the 
RF field efficiently drives single-quantum transitions but creates very little long-lived singlet component. 



In the low-RF-power regime, the long-lived coherence pllc is well-described by Equation 5 as a population difference 
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between the two central bare spin-pair eigenstates. The conventional two-spin dipole-dipole relaxation model of 
Solomon [551 predicts that in most cases Tuc — 3T X (see supplement SI). 



For intermediate RF spin- locking power (y n w Av), a more complex analysis is required. For an arbitrary RF power, 
the initial state psr can be represented as 

PsT=\n)(it\ + \it)(n\ 

+\{\<P+){<P + \ + \tf>+)(4>-\ + \4>-)(<t>+\ + Itf-W-I) 

-sin 2 0|0 s )(0 S |, (16) 



where the four spin-locked eigenstates are given by Eq. (7)-(10). 

At moderate RF powers (y n > Av), psr is still mainly composed of the population \4>s){4's\ an d mixed triplet 
states. However, the eigenstate \<f>s) no longer consists solely of the singlet \So). It also contains a triplet component 
cos0(|TL) + |T + ))/v2, which interacts with | <^>o) via a double-quantum transition, with relaxation rate scaling as 
cos 2 9. The triplet component also interacts with \4>+) and |</>_) via single-quantum transitions, with relaxation rate 
scaling as cos 2 9; and via double quantum transitions, with relaxation rate scaling as cos 2 9 sin 2 9. 

The above scaling of the relaxation of psr suggests a model for the measured singlet lifetime as a function of RF 
spin-locking power: 

1 1 2 a , 1 



.measured T x 

1 1 



cos" 11 { — (17) 

(18) 



T x l + (2v n /Av) 2 T s 

where T x is the lifetime at low RF power and Ts is the maximum singlet lifetime, typically achieved at high RF 
power. Significantly, this model predicts that for typical maximum singlet lifetimes, the measured singlet lifetime 
reaches 95% of its maximum value when the nutation rate v n is approximately 5 Av. 

Relaxation of the long-lived coherence pllc can be modeled using a similar analysis. In terms of the spin-locked 
eigenstates, we have: 

p L Lc=\n){n\-\m^\ 

= costf(|0 + )(0 + |-|0-X<M) 

+ ^(|0s)(^ + | + |0 + >^s| + \<t>s){<l>-\ + \<f>-)(4>s\)- (19) 



At low RF spin-locking powers, the long-lived coherence is mainly composed of |0 + )(<^ + | — |</>_)((^_|, and these two 
eigenstates interact with one another via a zero-quantum transition. However, as the RF power is increased, these 
states begin to mix with \T + ) and |T_), which opens up double-quantum transitions with relaxation rates scaling 
as sin 4 #. A double-quantum transition with |0q) also becomes available, with relaxation rate scaling as sin 2 9. The 
latter relaxation rate increases more quickly with RF power and dominates at small 9. 
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This above scaling suggests a simple model for the pllc relaxation rate: 



1 • 2 «> 1 (20) 



= — sm _ 

Tllc .measured -*-y Tllc 

1 (2^/A^) 2 1 

T y l + {2v n /^ v f + Tllc 1 ' 

where 1 /Tllc is the relaxation rate at zero RF power and 1/T y is the additional relaxation rate due to the applied 
RF power. 

We find that our model for the measured singlet lifetime agrees well with the detailed treatment of Pileio and 
Levitt (see eq. 43 in |18j). which contains terms up to eighth power in cos 9. Our model includes only lowest-order 
terms, but satisfactorily describes the measured relationship between singlet lifetime and RF power, as described 
below. The two models deviate most at low RF powers (y n < 0.5A^), where higher-order terms in cos 9 make larger 
contributions. See supplement S2 for a comparison. 



3. Experimental Results 



We performed NMR studies at 4.7 T of proton pair singlet states in a number of small organic molecules using a wide 
range of RF spin-locking powers. We chose citric acid and p-hydroxybenzoic acid, as Pileio et al. had previously 
studied these using high RF power [5]. Aditionally, we studied aspartic acid, trans-l,4-cyclohexanediol, and glycerol 
formal as examples of molecules with a range of structures. 



Bi (pT) Bj (uT) Bj 

5 10 15 02468 10 02468 10 

5 i , , , 1 1 , 1 , , 1 1 , , , 1 ^ 15 -) 2.2 - J 




Figure 3: Measurements of the lifetimes of the singlet state, long-lived coherence (LLC), and a mixture of the two as a function of the 
effective RF spin-locking field B± for proton pairs in small organic molecules. Also shown are fits to models for the singlet and LLC 
lifetimes, as described in the main text. (A) Citric acid, T x = 500 ms and T y = 600 ms; (B) p-hydroxybenzoic acid T x = 2 s and T y = 7 s; 
(C) glycerol formal, T x = 250 ms and T y = 1 s; RF power is quantified by the induced nutation frequency about the Bi field. Molecular 
structures are shown, protons of the singlet(s) are circled, and values for chemical shifts (At/) and spin-lattice relaxation times (Ti) are 
indicated. 

Our experimental protocol (Fig. 2) initialized proton pairs into one of three different states: the singlet-enhanced 
superposition state psT: which contained predominantly singlet population after rapid initial decay of the triplet 
component; a long-lived coherence between singlet and triplet, pllc'i or a mixture of the two. In all molecules we 
measured the singlet (pst) lifetime to increase with the applied RF power, reaching a plateau at the maximum singlet 
lifetime, Tg, when v n > 5Aia In contrast, we found that the lifetime of the long-lived coherence decreases from its 
maximum value of Tllc with the application of RF power. Both of these cases are well-modeled by equations 18 and 
21 above. When we created a mixture of psr and pllc-, the measured lifetime was that of the state with the longest 
lifetime at a given RF power, although the measured amplitude was lower as the contribution from the faster-relaxing 
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state was quickly lost. We individually fit the two regimes of the mixed-state lifetime-vs.-RF-power measurements 
with the corresponding models for pg? an d Pllc-, which provided a good characterization of the system's behavior, 
as shown in Fig. 3A-C. Note that for each molecule studied, we found that the shortest mixed-state lifetime occurs 
near v n 0.5A^. Results for maximum singlet and LLC lifetimes (T$ and Tllc) are summarized in Table 1. 



Table 1: Measured values of spin-lattice, singlet, and long-lived coherence (LLC) relaxation times. 

Molecule T x {s) T s {s) Enhancement (T s /Ti) T LLC {s) Enhancement {T LLC /T{) 

citric acid 0.58 ± 0.03 4.5 ± 0.3 7.8 ± 0.7 1.5 ±0.1 2.6 ± 0.2 

p-hydroxybenzoic acid D 2 2.9 ±0.1 16 ± 2 6.2 ± 0.8 7.3 ± 0.7 2.5 ± 0.3 

p-hydroxybenzoic acid H 2 2.3 ±0.1 5.8 ± 0.2 2.5 ±0.1 3.9 ±0.1 1.7 ±0.1 

aspartic acid 0.83 ± 0.03 7.48 ± 0.3 9.0 ± 0.5 2.3 ± 0.4 2.8 ± 0.5 

trans-l,4-cyclohexanediol 1.34 ±0.02 3.9 ± 0.9 2.9 ± 0.7 2.90 ±0.01 2.16 ±0.03 

glycerol formal 0.68 ± 0.01 1.91 ±0.03 2.81 ± 0.06 1.37 ± 0.02 2.01 ± 0.04 



We also investigated possible mechanisms for proton-pair singlet relaxation. First, we compared our measurements 
with those at other magnetic fields to probe the importance of chemical-shift anisotropy (CSA) relaxation (see Table 
1). Our result of T$ = 4.5 s and 7.8-fold lifetime enhancement over T\ for citric acid at 4.7 T are consistent with a 
previous measurement of 4.81 s and 7.6 T\ at 9.4 Tesla [5]. Due to hardware limitations, there was insufficient RF 
power to reach the maximum singlet lifetime for p-hydroxybenzoic acid. Nevertheless, a fit to data at finite RF power 
gave T$ = 16 s and 5.5 Ti, which is consistent with previous measurements [5]. Since these singlet-state lifetimes 
do not significantly depend on magnetic field strength, we conclude that CSA is not a primary mechanism of singlet 
relaxation. To distinguish between the spin-rotation and magnetic quadrupole relaxation mechanisms, we performed 
singlet lifetime measurements at several sample temperatures. We found that the singlet lifetime increases linearly 
with temperature, which identifies the magnetic quadrupole relaxation mechanism as dominant in such proton-pair 
singlet state molecules (see supplement S3). 



Note that most previous proton-pair singlet measurements were conducted using deuterated solvents, which should 
result in weaker singlet-solvent interactions and larger enhancements of singlet state lifetime. To test whether such 
lifetime enhancement changed in a normally protonated solvent, we studied p-hydroxybenzoic acid in both D2O and 
H2O. We found that the enhancement of both the singlet and LLC lifetimes were significantly lower in H2O (see 
Table 1). The enhancement is likely higher in D2O due to the substitution of deuterium for the phenolic proton as 
well as reduced dipolar interactions with nearby solvent protons. 



4. Discussion 



The above experimental results and associated modeling establish an operational spin-locking condition is° pt rj 5A^ 
to realize maximum singlet lifetime with minimal RF power. In the context of this operational condition, we can 
reassess the past work by Levitt and colleagues using high-power RF spin-locking fields [?J [TS] . As shown in Table 
2, most of the previous experiments employed v n » 5Ai/: i.e., they used much higher RF power than was needed to 
achieve a long singlet-state lifetime. For example, for citric acid Pileio et al. [5] used v n = 3.5 kHz, whereas /S.v = 72 
Hz at 9.4 T, which is an order-of-magnitude higher spin-locking field than necessary. 

Furthermore, we note that similarly low RF powers will be required for practical in vivo singlet-state creation in 
a wide variety of molecules using clinical MRI scanners, where the static magnetic field is commonly between 1.2 
and 7 T; see example values for v° pt at 1.5 T given in Table 2. For example, at 4.7 T glycerol formal's protons 
have a frequency difference Av « 16 Hz; hence v° pt 80 Hz is sufficient to achieve significant singlet-state lifetime 
enhancement. For common biomolecules such as citric acid and aspartic acid, v° pt < 100 Hz at 1.5 T, which is well 
within the spin-locking regime commonly used in clinical MRI [S7J 12H] ■ The spin- locking times and strengths used 
in [37J |2E] imply that a 60 Hz spin- lock could be safely applied for 3.5 s, and a 20 Hz spin- lock for 30 s. These 
timescales are of the same order as the singlet lifetimes we measured for typical small molecules, and thus should be 
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sufficient to conduct a variety of in vivo measurements using singlets. Alternatively, the long-lived coherence can be 
utilized without the need for any spin-locking if only moderate lifetime enhancements are required. 

In summary, our measurements and theoretical description show that for many molecules long-lived nuclear-spin 
singlet states and singlet /triplet coherences can be created using RF spin- locking powers that are more than two 
orders of magnitude lower than in previous studies; and that the effectiveness of the spin-locking can be accurately 
predicted from spectral parameters. These insights will be useful in the development of new applications for singlet 
states in vivo, where the RF specific-absorption rate (SAR) must be minimized. 



Table 2: Comparison of the chemical shift; optimal spin nutation frequency for RF spin-locking, RJ 5Af; and values of v n used in 
previous experiments. Also listed are values for Vn Pt for a clinical MRI scanner. 

Molecule Av (Bq field) is° pt v n used in previous experiments v° pt at 1.5 T 

citric acid [5] 72 Hz (9.4 T) 360 Hz 3500 Hz 57Efc 

p-hydroxybenzoic acid [5 445 Hz (11.75 T) 2224 Hz 3500 Hz 284 Hz 

aspartic acid [9] 100 Hz (11.75 T) 500 Hz 2500 Hz 64 Hz 



5. Materials and Methods 



Solutions of citric acid, aspartic acid, p-hydroxybenzoic acid, and 1,4-cyclohexancdiol were made in D2O, with the 
addition of sodium hydroxide where necessary for dissolution. Glycerol formal was analyzed neat. Concentrations 
and conditions can be found in Table 3. All reagents were purchased from Sigma- Aldrich . All samples were prepared 
in 10 mm diameter NMR sample tubes and bubbled with nitrogen gas for three minutes. Spectra were acquired on 
a 200 MHz Bruker spectrometer without spinning. 

Experiments shown in Fig. 2 were run on each compound using varying lengths for T4. Pulse sequence parameters can 
be found in Table 4. To remove any remaining triplet polarization, phase cycling was used in which the experiment 
was repeated with both the first and last 90° pulses along -x rather than x. Between 8 and 32 averages were used 
to provide sufficient signal-to-noise. The intensity of each peak was then measured and plotted against r 4 . The 
resulting data was fit with a single exponential time decay. Further details of the experimental pulse sequences are 
discussed in Supplement S4. Multiple datasets were collected using different RF power levels for spin-locking. The 
RF power was characterized by measuring the nutation rate induced by the RF B\ field, which was calibrated using 
a sequence of single-pulse experiments performed with increasing pulse length. T% relaxation rates were measured 
through conventional inversion-recovery experiments. 

For variable temperature experiments, the temperature was controlled by supplying hot air to the probehead. Blown 
air was heated with a Hotwatt cartridge heater controlled by an Omron temperature control box, and the temperature 
of the sample was monitored with an RTD in the probehead. 



Table 3: Sample preparations for the study of long-lived states. 

Molecule Concentration NaOH Concentration Solvent 



citric acid 0.26 M D 2 

p-hydroxybenzoic acid 0.29 M 0.50 M D 2 0, H 2 

aspartic acid 0.020 M 1.0 M D 2 

1,4-cyclohexanediol (cis/trans mixture) 0.41 M D 2 

glycerol formal neat neat 
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Table 4: Delays, in ms, for pulse sequences used in the experiments: ri — T2 — tz — r 4 — T5 

Molecule Singlet Coherence Mixture 

citric acid 12.5-3.7-6.0-r 4 -6.0 12.0-17.0-r 4 -15.5 12.0-12.0-r 4 -12.0 

p-hydroxybenzoic acid 29.3-31. 5-1. 14-r 4 -l. 14 2.0-3. 75-t 4 -1. 25 30.0-1. 25-r 4 -1.25 

aspartic acid 11.8-8.0-3.5-r 4 -4.0 - 7.4-7.0-t 4 -3.7 

1,4-cyclohexanediol (cis/trans mixture) - - 4.2-3.0-r 4 -2.0 

glycerol formal - - 31.0-15.0-r 4 -15.0 
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